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1 Field Of The Invention 

2 The invention relates generally to the field of computer graphics, and more specifically to 

3 generation of meshes representing surfaces of objects. The invention specifically provides a 

4 computer graphics system and computer-implemented method for generating a coarse-level mesh 

5 from a fme-level mesh, which may be used in multi-resolution surface editing. 

6 Background Of The Invention 

7 Generally, in computer graphics, obj ects are represented as surfaces, with the surfaces being 
J'pfS represented by meshes. A mesh consists of a set of vertices, or points in three-dimensional space, 
J5f9 interconnected by edges. The edges define polygonal faces, which may be in the form of triangles, 
fiO quadrilaterals, and so forth. In some computer graphic operations, it is desirable to generate a 
Ml representation of a surface at a finer resolution than a current representation. There are several 

2 popular methodologies for generating a representation of a surface at a finer resolution than a current 

Q3 representation, including a CatmuU-Clark surface subdivision methodology and a Loop surface 

:1| 4 subdivision methodology. In other operations, it can also be desirable to generate from a mesh 

^ 5 representing a surface at one level of resolution, a mesh at a coarser resolution. A methodology for 

lI 6 generating a representation of a surface at a coarser resolution than a current representation, and a 

17 methodology for generating such a surface representation, based on the Loop methodology, is 

18 described in D. Zorin, et al, "Interactive Multiresolution Mesh Editing, SIGGRAPH 1997, 

19 Computer Graphics Proceedings, pp. 259-268, particularly p. 262. 

20 The methodology described in Zorin relies on a method for smoothing triangular meshes that 

21 is described in G. Taubin, "A signal processing approach to fair surface design," SIGGRAPH 1995, 

22 Computer Graphics Conference Proceedings, pp. 351-358. Taubirfs smoothing methodology 

23 consists of two Unear filtering steps, the first being a low-pass filter step, followed by a high-pass 

24 filter step. The high-pass filter step is used to cancel out a mesh- shrinking effect that is induced by 

25 the low-pass filter step. The filtering operations described by Taubin have the problem that, as a 

26 consequence of the linearity of the filtering procedure, artifacts referred to as Gibbs ripples are 
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1 introduced in the smoothed mesh as a result of discontinuities in the input mesh. When used in 

2 generating a coarse level mesh, Taubin's filtering produces Gibbs ripples at the places on the surface 

3 where the fine level mesh has sharp variations. Although Taubin's smoothing methodology may be 

4 configured with the help of parameters, no filter can play the role of an inverse of the Loop 

5 subdivision rule. Accordingly, in a multi-resolution representation of a Loop subdivision surface, 

6 there will be redundant detail information. 

7 Summary Of The Invention 

8 The invention provides a new and improved system and method for generating a coarse-level 
...9 mesh from a fme-level mesh, which may be used in multi-resolution surface editing. 

CiO In brief summary, the invention provides an arrangement for generating a coarse ("j level 

f|l mesh representing a surface, from a finer (**j+l") level mesh surface representation. The 

li2 arrangement includes an indicator value generator and a coarse level mesh generator. The indicator 

013 value generator, for respective ones of the points in the finer level mesh surface representation, 

ri4 evaluates an indicator function, the value indicating whether a subdivision-inverse filter 

Hi 5 methodology or a least-squares optimization methodology is to be used to determme a position for 

ul6 the corresponding point in the coarse level mesh representation. The coarse level mesh generator 

y 7 determines, for each of the points that is to be provided in the coarse level mesh representation, a 

1 8 position in response to the position of the corresponding point in the finer level mesh representation, 

19 in accordance with the one of the subdivision-inverse filter methodology or least-squares 

20 optimization methodology indicated by the indicator value generated by the indicator value 

21 generator. 

22 Brief Description Of The Drawings 

23 This invention is pointed out with particularity in the appended claims. The above and 

24 further advantages of this invention may be better understood by referring to the following 

25 description taken in conjunction with the accompanying drawings, in which: 
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1 FIG. 1 depicts a computer graphics system including an arrangement for generating a coarse 

2 level mesh representation of a surface from a finer level mesh of a surface, constructed in accordance 

3 with the invention; 

4 FIG. 2 depicts a mesh representing a surface to which the Loop surface subdivision 

5 methodology has been applied; 

6 FIG. 3 depicts a mesh representing a surface to which the CatmuU-Clark surface subdivision 

7 methodology has been appUed; 

8 FIGS. 4 and 5 depict diagrams useful in understanding operations performed by the fme-to- 

9 coarse level mesh generating arrangement; 

i JO FIG. 6 depicts a flow chart describing operations performed by the fine-to-coarse level mesh 

?|l generating arrangement in connection with triangular meshes; and 

l|2 FIG, 7 depicts a flow chart describing operations performed by the fme-to-coarse level mesh 

ni3 generating arrangement in connection with quadrilateral meshes. 

i| 4 Detailed Description of an Illustrative Embodiment 

p5 FIG. 1 depicts a computer graphics system 10 including an arrangement for generating a 

16 coarse level mesh representation of a surface from a finer level mesh of a surface^ constructed in 

17 accordance with the invention. With reference to FIG. 1, the computer graphics system includes a 

18 processor module 1 1, one or more operator input devices 12 and one or more display devices 13. 

19 The display device(s) 13 will typically comprise a frame buffer, video display terminal or the like, 

20 which will display information in textual and/or graphical form on a display screen to the operator. 

21 The operator input devices 12 for a computer graphics system 10 will typically include a pen 14 

22 which is typically used in conjunction with a digitizing tablet 15, and a trackball or mouse device 

23 16. Generally, the pen 1 4 and digitizing tablet will be used by the operator in several modes. In one 

24 mode, the pen 14 and digitizing tablet are used to provide updated shading information to the 

25 computer graphics system. In other modes, the pen and digitizing tablet are used by the operator to 

26 input conventional computer graphics information, such as line drawing for, for example, surface 
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1 trimming and other information, to the computer graphics system 10, thereby to enable the system 

2 1 0 to perform conventional computer graphics operations. The trackball or mouse device 1 6 can be 

3 used to move a cursor or pointer over the screen to particular points in the image at which the 

4 operator can provide input with the pen and digitizing tablet. The computer graphics system 1 0 may 

5 also include a keyboard (not shown) which the operator can use to provide textual input to the 

6 system 10. 

7 The processor module 1 1 generally includes a processor, which may be in the form of one 

8 or more microprocessors, a main memory, and will generally include one a mass storage subsystem 

9 including one or more disk storage devices. The memory and disk storage devices will generally 

1 0 store data and programs (collectively, "information") to be processed by the processor, and will store 
fci processed data which has been generated by the processor. The processor module includes 

11 connections to the operator input device(s) 12 and the display device(s) 13, and will receive 
111 information input by the operator through the operator input device(s) 12, process the input 
1[4 information, store the processed information in the memory and/or mass storage subsystem. In 

addition, the processor module can provide video display information, which can form part of the 

1 6 information obtained from the memory and disk storage device as well as processed data generated 

if thereby, to the display device(s) for display to the operator. The processor module 1 1 may also 

1^ include connections (not shown) to hardcopy output devices such as printers for facilitating the 

p generation of hardcopy output, modems and/or network interfaces (also not shown) for connecting 

20 the system 10 to the public telephony system and/or in a computer network for facilitating the 

21 transfer of information, and the like. 

22 The invention provides an arrangement for generating a coarse-level mesh from a finer level 

23 mesh. In the following, the fine level mesh will continue to be identified by index "j+1," and the 

24 coarse level mesh will continued to be identified by index "j." Before describing the operations 

25 performed by the fine-to-coarse mesh generating arrangement in connection with generating a 

26 coarse level mesh from a current representation, it would be helpful to describe two surface 

27 subdivision methodologies, namely the aforementioned Loop surface subdivision methodology and 

28 a CatmuU-Clark surface subdivision methodology, which are used to generate a finer ("j+1) level 

29 mesh from a coarse ("j") level mesh. The invention provides an arrangement for essentially 
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performing, in the case of a triangular mesh, operations that generally provide the inverse of Loop's 
surface subdivision methodology, and, in the case of a quadrilateral mesh, operations that generally 
provide the inverse of the CatmuU-Clark surface subdivision methodology. 

Loop's surface subdivision methodology will be described in connection with FIG. 2. 
Generally, in Loop's surface subdivision methodology, each triangular face in the original mesh is 
split into a plurality of subface, the subfaces defining the fmer level mesh. The vertices of the finer 
level mesh are positioned using weighted averages of the positions of the vertices in the original 
mesh. More specifically, and with reference to FIG. 2, a mesh 20 at a fine "j" level is depicted that 
includes a vertex v^ 21(0) located at position d(q), and a plurality "K" of surrounding points v (1) 
through Vq(K) (generally, Vq(k)) 21(1) through 21(K) (generally identified by reference numeral 
21(k)) located at positions at respective positions c'(l) through c'(K), the points v^(k) 21(k) being 
those points in the mesh 20 that are connected to vertex v^ by edges. ("K," the number of points 
connected to vertex v^ will sometimes be referred to as the vertex's "valence.") Denoting the set of 
the indices of the points Vq(k) in the mesh that are connected to vertex v^ by N(q,j), a mesh at the 
next higher fineness level "j+1" corresponds to vertex v'^ 22(0) located at position c'^^(q) and a set 
of surrounding points v\(l) 22(1) at located at positions d*'(l) and connected to the vertex v'^ 22(0). 
The finer level mesh is constiiicted by providing the vertex v'^ 22(0) at a position c''''(q) that is 
determined by 

and points v'q(l) 22(1) at positions ^"^'(1) that are determined by 

^^^H/) = l\ciq)^c^{k)Y\[c^{k-\)^c^{k^\)\,k=l,...,K (2), 



where the weighting factor a(K) is given by 
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a(K) = - 



3 + 2 cos — 

V K) 

8 



(3). 



1 

2 Equations (1) and (2) are applied by selecting each point in the mesh at level "j" as a vertex, to 

3 provide points for the mesh at level ." Loop's surface subdivision methodology can be applied 

4 recursively to provide a mesh at any desired degree of fineness. It will be appreciated that equation 
\5 (2) will provide that the point v'q(l) 22(1) that is created when the methodology is applied to point 
7^ Vq 21(0) as the vertex will be m the same position when the methodology is applied to point Vq(k) 
yt 2 l(k) as the vertex. 

p| It will be appreciated that, the Loop surface subdivision methodology provides 

s 9 (i) one new point, at location ^^^(1) in the finer ("j+l") level mesh, for, and generally 

IP somewhere between, each pair of points in the coarse ("j") level mesh, and 

|l (ii) a relocated point, at location c'^^(q) in the finer level mesh for, and generally 

iL2 somewhat near, each point in the coarse ("j") level mesh, 

13 which are interconnected to form the finer ("j+1) level mesh. Accordingly, in performing the 

14 inverse, that is, processing a finer ("j+1") level mesh to generate a coarse ("j") level mesh, the 

15 arrangement processes the points in the finer ("j+1") level mesh to remove, for each pair of 

1 6 connected points in the finer ("j+1 ") level mesh, and determines a position for the remaining point, 

1 7 which will be present in the coarse ("j ") level mesh, based on the positions of surrounding points in 

18 the finer ("j+1") level mesh. The operations performed by the arrangement in determining the 

19 position of the point that is to remain in the coarse ("j") level mesh will be described below, 

20 The Catmull-Clark subdivision methodology produces smooth surfaces using a small number 

21 of neighboring vertices. The Catmull-Clark surface subdivision methodology will be described in 
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1 connection with FIGS. 3A through 3E. With reference to FIG. 3A, that FIG. depicts a mesh 30 

2 consisting of four quadrilaterals 31(1) through 3 1 (4) . Each quadrilateral is referred to as a face. The 

3 quadrilaterals are defined by points 3 2(0), which is common to all of the quadrilaterals 31(1) through 

4 31(4), and other points 32(1) through 32(8). Generally, the CatmuU-Clark surface subdivision 

5 methodology is performed in a series of iterations, including 

6 (i) a face point generation iteration, in which, for each face, a face point 33(1) through 33(4) 

7 is generated as the average of the points 32(p) defining the respective face (reference FIG. 3B); 

8 (ii) an edge point generation iteration in which, for each edge, a new edge point 34(1) 

9 through 34(4) is generated as the average of the midpoint of the original edge with the two new face 
JO points of the faces adjacent to the edge (reference FIG, 3C; the face points generated during the face 
if point generation iteration are shown as circles); 

1 2 (iii) a vertex point iteration in which a vertex point 3 5 is generated in relation to the positions 

if> of the new face points generated in iteration (i) for faces that are adj acent to the original vertex point, 

Jjl the positions of the midpoints of the original edges that are incident on the original vertex point, the 

45 position of the original vertex point, and the valence of the original vertex point (reference FIG. 3D; 

§6 the face and edge points generated during the face point and edge point generation iterations are 

^7 shown as circles); and 

p (iv) a mesh connection step (reference FIG. 3E) in which: 

1 9 (a) each new face point is connected to the edge points of the edges defining the original 

20 face; and 

21 (b) each new vertex point is connected to the new edge points of all original edges 

22 incident on the original vertex point. 

23 More specifically, in the CatmuU-Clark subdivision methodology, denoting, for a vertex v^ 

24 at location d(q), 

25 (i) the set of indices of neighboring points that are connected thereto (for example, for 

26 point 32(0), point 32(2) and 32(8)) by N,(q,j), and 
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(ii) the set of indices of second-order neighboring points that are opposite to vertex ^(q) with 
respect to the level "j" faces (for example, for point 32(0), point 32(1)) that are incident with vertex 
d(q)byN^q,j) 

a mesh at the next higher fineness level "j+1 " is constructed as follows. In the face point generation 
iteration, the face points are generated and located at positions c'"^^(mi) determined as follows 



In the edge point generation iteration, the edge points are generated and located at positions ^"^'(1,) 
determined as follows: 

G NXqJ^ l), A:, G NXqj),m,_^,m, g N^{q,j+ l) 



In the vertex point generation iteration, the new vertex points are generated and located at positions 
as follows: 



In terms of only points that is substituting for the last term in equation (6), 



(7). 
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2 Generally, for quadrilateral meshes, the arrangement makes use of the Catmull-Clark 

3 methodology, except for the case of K=3 in equation (6). In that case, the arrangement makes use 

4 of 

3 

6 where y = — . Equation (8) would correspond to the Catmull-Clark methodology (equation (6)) 
38 

Si 3 1 

^ with / = -^-^ , except that, for K=3, / - instead ^^^2' t^™^ ^^^y points c' (compare 
^ equation (7)) 

U ^'"i^) = (l-7/)c^(^)+6r^ Zc^(/) + 7-| Tc^im) (9). 



10 As is apparent from the above discussion, the Catmull-Clark surface subdivision 

1 1 methodology applied to a quadrilateral mesh at the coarse ("j") level will produce a quadrilateral 

12 mesh at the finer ("j+l") level having the same number of points as the mesh at the coarse ("j") level, 

13 but the points in the finer ("j+l") level may be at different positions than in the coarse ("j") level. 

14 Accordingly, in performing the inverse, that is, processing a finer ("j+l") level mesh to generate a 

15 coarse ("j") level mesh, the arrangement processes the points in the finer ("j+1") level mesh 

16 determine a position for each point in the finer ("j+l") level mesh based on the positions of 

17 surrounding points in the finer ("j+1") level mesh. The operations performed by the arrangement 

18 in determining the position of the point in the coarse ("j") level mesh will be described below. 

19 As noted above, the invention provides an arrangement for generating a coarse ("j) level 

20 mesh from a finer ("j+l") level mesh. In generating a coarse mesh from a fine mesh, an indicator 
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function is used locally in connection with each point that is being processed as a vertex to determine 
whether, around the respective point, the fine mesh can be derived by subdivision of a coarser mesh 
using a subdivision inverse filter. If so, the value of the indicator function, which will be described 
below, will be "zero" or close to "zero." On the other hand, if the fine mesh cannot be derived by 
subdivision of a coarser mesh in the region around the respective vertex using a subdivision inverse 
filter, the value of the indicator function will not be close to "zero," and in that case the coarse mesh 
will be generated using an approximate solution of a least-squares problem. 

The indicator function is constructed for each point, selected as a vertex v'^, using the value 
of the Laplacian L(k,j+1) of the positions of the point and of the other points in its neighborhood in 
the finer (j+l) level mesh. For triangular meshes, in the one-dimensional case, for example, for a 
vertex v'^ on a boundary, crease line or the hke in the finer ("j+l ") level mesh, the Laplacian for the 
vertex is 



(10), 



where c'^^(k) is the position of the vertex v ^ in the finer ("j+l) level mesh for which the Laplacian 
is being generated, and c^"" ^ (k- 1 ) and d""^ (k+ 1 ) are the positions of the neighboring points v ^(k- 1 ) and 
v*^(k+l) in the fine (j+l) level mesh. In the two-dimensional case, that is, in the case of a vertex v'^ 
other than on a boundary, crease line or the like, for a regular vertex, that is, for a vertex for which 
the valence "K" is equal to "six," the Laplacian is 

if 1 

L{kJ^\) = - Zc^^H/) -c''\k) (11), 



\leN{kj + l) 



where the factor "1/6" reflects the valence of the vertex. Also in the two-dimensional case, for an 
irregular vertex, that is, for a vertex for which the valence "K" is not equal to "six," which may be 
the case at, for example, edges of the mesh, the Laplacian is 
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L{k,j+\) = p 



^ leN{kj+\) 



(12), 



1 

2 

3 



where 



3+8 a{K) 
^ ' ~ 3{-5 +Sa(K)) 



(13), 



m 



where a(K) is defined in equation (3) above. 

For quadrilateral meshes, in the one-dimensional case, that is, for a vertex v'^ on a boundary, 
crease line, or the like in the finer ("j+1") level mesh, the Laplacian L(k,j+1) for the vertex v^ is 



(14) 



To 
11 

12 
13 
14 

15 



where (^"^'(k) is the position of the vertex v'^ for which the Laplacian is being generated in the finer 
("j+1") level mesh, and ^""(k-l) and ^^'(k+l) are the positions of the neighboring points in the finer 
(j+1) level mesh. In the two-dimensional case, that is, in the case of a vertex v'q other than on a 
boundary, crease line or the hke, two Laplacians are used, namely a Laplacian Lg(kj+1) using the 
position d*'(k) of the vertex v'^ and the positions of the points whose indices are in the set Ne(qj), 
and a Laplacian Lf(kJ+l) using the position c''"(k) and the positions of the vertices whose indices 
are in the set Nf(q,j), as follows 



_1_ 
~K 



,7+1 



(15) 



16 



and 
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'k 



( 



- c 



7+1 



(16). 



\lsN^{k,j+i) 



As will be described below, the arrangement makes use of values of both Laplacians Lg(k,j+1) and 
Lj(kJ+l), and so both values will be generated for each vertex. 

As noted above, the arrangement makes use of an indicator function. For triangular meshes, 
the indicator function i(q) for a point as vertex v'^ is given by, in the one-dimensional case 

i{q) = L{q)-{L{q-\)+L{q^l))+^{L{q-2)+L{q+2)) (17) 
and, for the two-dimensional case 

i(q) = -6L{q)+2 S - 2 2^(0+ 1 L{l) ^^^X 

leN(\,0){q) leN{2,0){q) leN(2,l)(g) 

where the indices for the sums in equation (18) are in the neighborhood of the vertex for which the 
indicator function is being evaluated are as depicted in FIG. 4. 

In generating the values of the indicator function for quadrilateral meshes, the arrangement 
makes use of a linear combination L-^Jk,]-^!) of the Laplacians L^Ckj+l) and Lf(kj+1) described 
above (equations (15) and (16), respectively), namely 

LjkJ+\) = 2^LXk,j+l)+^L,{kJ+l) (19) 

for vertices for which valence "K" is not equal to "three." For vertices for which the valence "K" 
does equal "three," Lij^^j(kJ+l) takes the form 
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A„.(^.y+l) = - — LXkJ^\)^2Lf[k,j^\) (20). 

2 For regular vertices, that is, for vertices for which the valence "K" equals "four/' equation (19) 

3 becomes 

4 LZ = -2LXkJ+l)+ L,{k,j+l) (21). 

5 With this background, the indicator function i(k) for quadrilateral meshes is given by (omitting the 

6 index "j+1 since the indicator function will always be evaluated for vertices in the fine (j+1) level 
Ol mesh) 

if i{k) - jljk)- lLjl),\ iLjl), iLjl) 

S 111 

^ ^ /€A^(3,l)a) ^ /eiV(3,2)U) ^ /eA^(4,2)(^) 

UI9 where the vertex neighborhoods for the respective sums are as shown in FIG. 5. 

i 0 As noted above, the fme-to-coarse mesh generating arrangement generates a coarse ("j ") level 

1 1 mesh from a finer ("j+l") level mesh. For a triangular mesh, the fine-to-coarse mesh generating 

12 arrangement determines, for each vertex v'^ located at position c'^^(k) in the finer ("j+l") level mesh 

13 for which a vertex v^ is to appear in the coarse ("j") level mesh, a position d(k) for the vertex v^ in 

14 accordance with 

1^ c'{k) = c''\k)-\- A L(k,J+\) (23), 

16 where ?i is a parameter whose value are determined as described below and L(k,j+1) is the value of 

17 the Laplacian generated for the vertex v^ in the finer (j+1) level mesh. As noted above, not all of 
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1 the vertices in the finer ("j+1 ") level mesh will have corresponding vertices in the coarse ('7') level 

2 mesh. If the value of the Laplacian L(kj+1) generated for a vertex v'^ in the finer ("j+l") level mesh 

3 is relatively small, that is, if 

|l(^,7>i)|<^ (24) 

5 for a selected small 8, then for equation (23), the position d(k) of the corresponding vertex v^ in the 

6 coarse ("j") level mesh can be taken as the position d''\k) of the vertex v'^. In that case, neither the 

7 parameter X, nor the indicator function i(q) (equations (17) and (18)), need to be generated. 

8 If the value of the Laplacian L(k,j+1) is not small, that is, if the condition in equation (24) 
ifp does not hold, the manner in which parameter X for use in equation (23) is to be determined will 
Sp depend on the value of the indicator fimction, generated as described above in connection with 
Uk equations (17) and (18) and whether the vertex is on a boundary, crease line or the like, a regular 
m vertex (that is, a vertex v'^ whose valence "K" is "six") or an irregular vertex (that is, a vertex v ^ 
% whose valence "K" is other than "six"). If the value of the indicator function is zero, or close to zero, 
G:4 for the vertex v'q, a subdivision-inverse filter is used in which 

J|S (i) if the vertex is on a boundary, crease line or the like 

S6 A = -I (25), 

17 (ii) if the vertex v'^ is regular, that is, if its valence "K" is equal to "six," and is not on a 

18 boundary, crease line or the like 

19 3 

i = -- (26), and 

2 

20 (iii) if the vertex v'^ is irregular, that is, if its valence "K" is not equal to "six," and is not on 

21 a boundary, crease line or the like 
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X = 



-5+ Sa{K) 



(27), 



where a(K) is as defined above in connection with equation (3). 

On the other hand, if the value of the indicator function i(q) is not zero or close to zero, a 
least-squares optimization is used in which 

(i) if the vertex is on a boundary, crease line or the like 



Z = 



L{k) 



(28), 



12 



where - - and = - 



23 
49 



(ii) if the vertex v'^ is regular, that is, if its valence "K" is equal to "six," and the vertex is not 
on a boundary, crease line or the like 



A = 



1 



Lik) 



" l^N{kJ+\) 



(29), 



61 14403 1 

where K^^ - and = ; it will be appreciated that the factor " — in equation 

^ 5720 ^ 5720 6 

(29) reflects the value of the valence "K," namely, "six." 

(iii) if the vertex v ^ is irregular, that is, if its valence "K" is not equal to "six," and the vertex 
is not on the boundary, crease line or the like 
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L{k) 



(30), 



where 



j^irreg 
^0 



2 (5- Sa{K)) (l4647 K - 391848 a(K) + 391848 
715(3+ Sa{k)){256+ 4\K- 512 a{k) + 256 a{kY) 



(31) 



^^4 and 

16(-5531i<:- 24521 a{K) + 24521 a{KY) 

rii z>r" = — ^ (32), 

' 715(256+ 4lK-5\2a{K) + 25a{Ky) 

life where, for both equations (3 1) and (32), *'K" corresponds to the valence of the irregular vertex, and 

Jriz a(K) is as defined in equation (3). 

h-^ For a quadrilateral mesh, for each vertex y\ located at position d''\k) in the finer G+1) level 

9 mesh, a vertex v^ is provided in the coarser (j) level mesh at a position d(k) that is determined by 



1 1 where and }^ are parameters whose values are determined as described below, and Le(k,j+1) and 

12 Lf(k,j+1) are determined as described above in connection with equations (14) through (16). In the 

1 3 one-dimensional case, that is, for vertices on a boundary, crease line or the like, since only Laplacian 

14 is used, only one parameter need be generated for use in equation (33), On the other hand for 

1 5 vertices that are not on a boundary, crease line, or the like, since both Laplacians and Lf are used, 

1 6 both parameters and are generated for use in equation (33). As with the triangular mesh, if the 
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values of L,(kj+1) and Lf(k,j+1) are relatively small, that is, if for a vertex v'^ in the fine ("j+1") 
level mesh 



max 



(\LXkJ+ll\L^{kJ+lfj<s 



(34) 



for a selected small e, then for equation (34), the position c'(k) of vertex v^ in the coarse ("j") level 
mesh can be taken as the position c^^'(k) of vertex v'^. In that case, neither the parameters and X^, 
nor the indicator function i(k,j+l) (equation (22)), need to be generated. 

If the condition in equation (34) does not hold, the values of parameter 1^ and for use in 
equation (33) will depend on the value of the indicator function. For a vertex v'^ that falls on a 
boundary, crease line or the like, the value of the parameter will be generated in the same maimer 
as in a triangular mesh (reference equations (25) and (28) above). On the other hand, for a vertex 
v'q that does not fall on a boundary, crease line or the like, if the value of the indicator function 
i(k,j+l) generated therefor is near zero, that is, if 



a subdivision-inverse filter is used in which, if the valence "K" of the vertex not equal to "three," 



e • max 



(|4(^,7 + 4|l^(^,7+i 



I 



(35) 



K- 3 



(36), 




K- 3 



and, if the valence "K" for the vertex is equal to "three,' 




(37). 
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On the other hand, if, for a vertex v'^, the value of the indicator function i(k,j+l) is not near zero, that 
is, if the condition in inequaUty (35) is not satisfied, 



1 



(38), 



and 



(39), 



which reflect a least squares optimization. If vertex v'^ is regular, that is, if the vertex's valence "K" 
is "four," 





9946871 




" 4862025 




1024 


^1 - 






~ 405 




1644032 




972405 




1338874 


^21 


' 972405 



(40). 



On the other hand, if the vertex v' is irregular, that is, if the vertex's valence "K" is other than "four," 
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162307143936- 92746939392^- 8924282387 

■LCC _ ^ 

~ 4862025 (l2544 - 14336 K + 4096 + 901 K') 
1024 (2793728 - 1596416 K - 244001 K^) 
" " 99225 (12544 - 14336 K + 4096 + 901 K^) 

5 12 (-1 13305472 + 64745984i^+ 17391149 ii:') 
4862025 (l2544 - 14336^ + 4096 + 901 K^) 
4(8660934688- 4949105536^- 1876158821^^) 
4862025(12544- 14336^+ 4096^' + 901 



(41). 



The following are approximations for equations (40) that are used in one embodiment of the 
arrangement: 



bll = - 2.04583 
Z>ff = - 2.5284 
b^, = 1.69069 
Z>2'[ = - 1.37687 

and the following are approximations for equations (41) that are used m the same embodiment: 



-19- 



MENT-060 



1623.07- 927.469 i:- 89.2428 J^' 
609.892 - 697.02 K + 199.149 + 43.8068 

28.6078 - 16.3473 K - 2.49857 
12.4468- 14.2249^+ 4.06426 A:' + 0.894017 

(43). 

-580.124 + 331.499 ^ + 89.0427 
609.892 - 697.02 K + 199.149 i?:^ + 43.8068 

346.437 - 197.964 i^- 75.0464 
609.892 - 697.02 K + 199.149 + 43.8068 

1 

C:i With this background, the operations performed by the fine-to-coarse mesh generating 

11 arrangement in generating a coarse ("j") level mesh from a finer ('^j+l ") level mesh will be described 
i!l in connection with the flow charts in FIGS . 6 and 7 . The flow chart in FIG. 6 depicts the operations 
i[B performed in connection with a triangular mesh and the flow chart in FIG. 7 depicts the operations 
ff| performed in connection with a quadrangular mesh. With reference initially to FIG. 6, the 
:7 arrangement will first perform a series of steps to generate the Laplacian (equations (10) through 
m (12)) for each of the points in the triangular finer ("j+1 ") level mesh. Initially, the arrangement will 
4^ select a point in the finer (j+1) level mesh as vertex (step 100) and determine whether the point is 
fib regular, that is, whether its valence "K" is six, and does not fall on a boundary, crease hne or the hke 
'l 1 (step 101). If the arrangement makes a positive determination in step 101, it proceeds to step 102 

12 in which it generates the Laplacian for the point selected in step 100 in accordance with equation 

13 (11). 

14 On the other hand, if the arrangement makes a negative determination in step 101, it will 

1 5 determine whether the point is irregular, that is, whether its valence "K" is less than six and does not 

16 fall on a boundary, crease line or the like (step 103). If the arrangement makes a positive 

17 determination in step 103, it proceeds to step 104 in which it generates the Laplacian for the point 

18 selected in step 100 in accordance with equation (12). 



^10 



^1 - 



^20 



b: 



21 
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1 Finally, if the arrangement makes a negative determination in step 103, which will be the 

2 case if the point falls on a boundary, crease line or the like, it proceeds to step 105 in which it 

3 generates the Laplacian for the point selected in step 100 in accordance with equation (10). 

4 Following step 102, 104 or 105, the arrangement will determine whether all of the points 

5 have been selected in step 1 00 (step 1 06). If the arrangement makes a negative determination in step 

6 106, it will not have generated Laplacians for all of the points in the finer ("j+1 ") level mesh, and, 

7 in that case, it will return to step 100 to select another point and perform respective steps 101 

8 through 1 05 in connection therewith. The arrangement performs step 1 00 and respective steps 1 0 1 

9 through 1 05 through a plurality of iterations, in each iteration selecting and generating the Laplacian 

1 0 for each of the points in the finer ("j+1 ") level mesh. 

11 After all of the points have been selected and Laplacians generated therefor, the arrangement 
II will make a positive determination in step 106, at which point it will sequence to a series of steps 
% in which the arrangement determines the positions c'(q) forpoints in the finer ("j+1") levelmesh that 
m are to be retained in the coarse ("j") level mesh, hi those operations, the arrangement will select a 
15 point in the finer ("j+1 ") level mesh (step 1 10) and determine whether the pomt is to be retained in 
% the coarse ("j") level mesh (step 1 1 1). If the arrangement makes a negative determination in step 
av 1 1 1 , the point is not to be retained in the coarse ("j ") level mesh. In that case, the arrangement will 

sequence to a step 1 12 in which it determines whether it has selected all of the points in the finer 

to ("j+1") level mesh. If the arrmgement makes a negative determination in step 1 12, it will return to 

20 step 11 0 to select another point in the finer ("j+1 ") level mesh. On the other hand, if the arrangement 

21 makes a positive determination in step 1 1 2, which will be the case if it has processed all of the points 

22 in the finer ("j+1 ") level mesh, it will exit (step 1 13). The arrangement can make use of any of a 

23 number of methodologies in determining whether a point in the finer ("j+l") level mesh will be 

24 retained in the coarse ("j ") level mesh. For example, given a suitable point indexing methodology, 

25 as will be apparent to those skilled in the art, the arrangement can determine whether a point in the 

26 finer ("j+1") level mesh will be retained in the coarse ("j") level mesh. 

27 Returning to step 1 1 1 , if the arrangement makes a positive determination in that step, which 

28 will be the case if the selected point in the finer ("j+1 ") level mesh will be retained, the arrangement 

29 sequences to step 1 15, in which it determines whether the value of the Laplacian for the point 
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1 selected in step 1 1 0 is less than "e" (reference inequality (24)). If the arrangement makes a positive 

2 determination in step 11 5, it determines the position c'(q) of the point in the coarse ("j ") level mesh 

3 as the position c'^'(q) of the corresponding point in the finer ("j+1") level mesh (step 116). 

4 Following step 1 1 6, the arrangement will sequence to step 1 1 2 to determine whether all of the points 

5 in the finer ("j+1 ") level mesh have been processed, and, if not, return to step 1 10 to select another 

6 point, and, if so, exit. 

7 Returning to step 115, if the arrangement makes a negative determination in that step, it will 

8 sequence to step 1 17, in which it will determine whether the point is regular, that is, whether its 

9 valence "K" is six, and does not fall on a boundary, crease line or the like. If the arrangement makes 

10 a positive determination in that step, it will generate a value for the indicator Amotion (reference 
l[| equation (1 8)) (step 118). The arrangement will then determine whether the value of the indicator 
m function is zero or close to zero (step 119). If the arrangement makes a positive determination in 
P step 1 1 9, it will provide a value for the parameter 1 as described above in connection with equation 
|l (26) (step 120). On the other hand, if the arrangement makes a negative determmation in step 1 1 9, 

11 it will provide a value for the parameter X as described above in connection with equation (29), 
16 provided that the value of the parameter I does not fall below the value specified in equation (26) 
i]l7 (step 121). 

Returning to step 1 1 7, if the arrangement makes a negative determination in that step, it will 

im sequence to step 122 to determine whether the point is irregular, that is, whether its valence "K" is 

20 not equal to six, and the point does not fall on a bovmdary, crease line or the like. If the arrangement 

21 makes a positive determination in that step, it will generate a value for the indicator function 

22 (reference equation (1 8)) (step 123). The arrangement will then determine whether the value of the 

23 indicator function is zero or close to zero (step 124). If the arrangement makes a positive 

24 determination in step 124, it will provide a value for the parameter 1 as described above in 

25 connection with equation (27) (step 125). On the other hand, if the arrangement makes a negative 

26 determination in step 124, it will provide a value for the parameter X as described above in 

27 connection with equation (30), provided that the value of the parameter 1 does not fall below the 

28 value specified in equation (26) (step 126). 
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1 Returning to step 1 22, if the arrangement makes a negative determination in that step, it will 

2 sequence to step 127 to determine whether the point falls on a boimdary, crease line or the like. If 

3 the arrangement makes a positive determination in that step, it will generate a value for the indicator 

4 function (reference equation (17)) (step 128). The arrangement will then determine whether the 

5 value of the indicator fimction is zero or close to zero (step 129). If the arrangement makes a 

6 positive determination in step 1 29, it will provide a value for the parameter X as described above in 

7 connection with equation (25) (step 130). On the other hand, if the arrangement makes a negative 

8 determination in step 129, it will provide a value for the parameter X as described above in 

9 connection with equation (28), provided that the value of the parameter X does not fall below the 

1 0 value specified in equation (25) (step 131). 

0 Following step 120, 121, 125, 126, 130, or 131 , the arrangement will sequence to step 132, 

11 in which it determines whether the value of the parameter X generated in the respective step was 
ll greater than "one," and, if so, limit the maximum value of parameter X to "one" (step 133), 
€1 Thereafter, the arrangement uses equation (23) to determine the position d^\q) of the point in the 
il coarse ("j") level mesh (step 134). 

ClB Following step 1 34, the arrangement the arrangement will sequence to step 11 2 to determine 

If whether all of the points in the finer ("j+l") level mesh have been processed, and, if not, return to 

'a| step 1 10 to select another point, and, if so, exit. 

^19 The flow chart in FIG. 7 depicts the operations performed in connection with a quadrangular 

20 mesh. With reference to FIG. 7, the arrangement will first perform a series of steps to generate the 

2 1 Laplacians (equations (14) through ( 1 6)) for each of the points in the quadrilateral finer ("j+l ") level 

22 mesh. Initially, the arrangement will select a point in the finer (j+l) level mesh as vertex (step 1 50) 

23 and determine whether the point does not fall on a boundary, crease line or the Uke (step 151), If 

24 the arrangement makes a positive determination in step 151, it proceeds to step 152 in which it 

25 generates the Laplacian for the point selected in step 1 50 in accordance with equations (1 5) and (1 6). 

26 On the other hand, if the arrangement makes a negative determination in step 151, which will 

27 be the case if the point falls on aboundary, crease line or the like, it proceeds to step 153 in which 

28 it generates the Laplacian for the point selected in step 150 in accordance with equation (14). 
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1 Following step 152 or 153, the arrangement will determine whether all of the points have 

2 been selected in step 150 (step 1 54). If the arrangement makes a negative determination in step 1 54, 

3 it will not have generated Laplacians for all of the points in the finer ("j+1 ") level mesh, and, in that 

4 case, it will return to step 1 50 to select another point and perform respective steps 1 5 1 through 1 53 

5 in connection therewith. The arrangement performs step 150 and respective steps 151 through 153 

6 through a plurahty of iterations, in each iteration selecting and generating the Laplacian for each of 

7 the points in the finer ("j+1 ") level mesh. 

8 After all of the points have been selected and Laplacians generated therefor, the arrangement 

9 will make a positive determination in step 1 54, at which point it will sequence to a series of steps 

1 0 in which the arrangement determines the positions c'(q) for respective points in the coarse ("j ") level 
G mesh using the positions d^\q) of the points in the finer ("j+1") level mesh and the values of 

11 parameters to be generated and the Laplacians generated in steps 150 through 153. Initially, the 
II arrangement will select a point in the finer ("j+1 ") level mesh (step 1 60) and determine whether the 
iii values of both Laplacians L,(kj+1) and Lf(k,j+1) generated for the point selected in step 1 60 are less 
ll than "£" (reference inequahty (34)) (step 161). If the arrangement makes a positive determination 
IP in step 161, it determines the position c'(q) of the point in the coarse ("j") level mesh as the position 
m c'^'(q) of the same point in the finer ("j+1") level mesh (step 162). FoUowmg step 162, the 
ft arrangement will sequence to step 163 to determine whether all of the points in the finer ("j+1") level 
is mesh have been processed, and, if not, return to step 1 60 to select another point. On the other hand, 
2^0 if the arrangement makes apositive determination in step 163, it will exit (step 164). 

2 1 Returning to step 1 6 1 , if the arrangement makes a negative determination in that step, it will 

22 sequence to step 1 65, in which it determines whether the determine whether the point is regular, that 

23 is, whether its valence "K" is four, and does not fall on a boundary, crease line or the hke. If the 

24 arrangement makes a positive determination in step 165, it will generate a value for the indicator 

25 function (reference equation (22)) (step 166) and determine whether the value of the indicator 

26 function is near zero (reference inequality (35)) (step 167). If the arrangement makes a positive 

27 determination in step 1 67, it will generate values for \ and as described above in connection with 

28 equations (36) (step 1 68). On the other hand, if the arrangement makes a negative determination in 

29 step 1 67, it will generate values for and as described above in connection with equations (3 8)- 
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1 (40), in the process ensuring that the values ofX^ and Xj are within predetermined Hmits (-4< A,; < 1/2, 

2 1/4<X2<1) (step 169). 

3 Returning to step 1 65, if the arrangement makes a negative determination in that step, it will 

4 sequence to step 170, in which it determines whether the point is irregular, that is, whether its 

5 valence "K" is other than four, and does not fall on a boundary, crease hne or the like. If the 

6 arrangement makes a positive determination in step 170, it will generate a value for the indicator 

7 function (reference equation (22)) (step 171) and determine whether the value of the indicator 

8 function is near zero (reference inequahty (35)) (step 172). If the arrangement makes a positive 

9 determination in step 172, it determine whether the valence "K" of the point is "three" (step 1 73). 

10 If the arrangement makes a positive determination in step 1 73 , it will generate values for and 

11 as described above in connection with equation (37) (step 174). On the other hand, if the 
II arrangement makes a negative determination in step 173, it will generate values for and as 
JjJ described above in connection with equation (36) (step 175). Returning to step 172, if the 
& arrangement makes a negative determination in step 172, which will be the case if the value for the 
II indicator function as generated in step 171 is not near zero, it will generate values for and as 
%p described above in connection with equations (38)-(40), m the process ensuring that the values of 

iii f 4 1 1 I] 

Bl7 Xi and X^ are within predetermined limits - ^ < X ^<—, — < A 2^ ^_ (step 176). 



1 8 Returning to step 170, if the arrangement makes a negative determination in that step, it will 

1 9 sequence to step 1 77 to determine whether the point falls on a boundary, crease line or the like. If 

20 the arrangement makes a positive determination in that step, it will generate a value for the indicator 

21 function (reference equation (17)) (step 178). The arrangement will then determine whether the 

22 value of the indicator fimction is zero or near zero (reference inequality (35) above) (step 179). If 

23 the arrangement makes a positive determination in step 1 79, it will provide a value for the parameter 

24 X, as described above in connection with equation (25) (step 180). On the other hand, if the 

25 arrangement makes a negative determination in step 179, it will provide a value for the parameter 

26 X as described above in connection with equation (28) (step 181). 
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1 Following step 168, 169, 174, 175, 180, or 181, the arrangement will sequence to step 1 82, 

2 in which it uses equation (33) and values for parameters and I2 generated in the previous steps to 

3 determine the position d^'Cq) of the point in the coarse ("j") level mesh. Following step 1 82, the 

4 arrangement will sequence to step 1 63 to determine whether all of the points in the finer ("j+ 1 ") level 

5 mesh have been processed, and, if not, return to step 160 to select another point, and, of so, sequence 

6 to step 164 to exit. 

7 The arrangement provides a number of advantages. In particular, the invention provides an 

8 arrangement for generating a coarse ("j ") level mesh representing a surface, from a finer ("j +1") level 

9 mesh surface representation. The arrangement makes use of an indicator function to provide an 

1 0 indication as to which of several methodologies is be used at respective ones of the points in the finer 
|:| ("j+1") level mesh, including a subdivision-inverse filter methodology and a least-squares 

optimization methodology. 

11 It will be appreciated that numerous changes and modifications may be made to the 
-ii arrangement as described herein. For example, it will be appreciated that, although the arrangement 
% has been described in coimection with meshes of triangular and quadrangular faces, it will be 
ttB appreciated that the arrangement may find utihty in connection with meshes having faces of different 
W polygonal structures. In addition, although equations for generating the indicator functions and 
^|:8 parameters may differ than the specific ones described herein. 

Y9 It will be appreciated that a system in accordance with the invention can be constructed in 

20 whole or in part from special purpose hardware or a general purpose computer system, or any 

2 1 combination thereof, any portion of which may be controlled by a suitable program. Any program 

22 may in whole or in part comprise part of or be stored on the system in a conventional manner, or it 

23 may in whole or in part be provided in to the system over a network or other mechanism for 

24 transferring information in a conventional manner. In addition, it will be appreciated that the system 
2 5 may be operated and/or otherwise controlled by means of information provided by an operator using 

26 operator input elements (not shown) which may be connected directly to the system or which may 

27 transfer the information to the system over a network or other mechanism for transferring 

28 information in a conventional manner. 
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The foregoing description has been Hmited to a specific embodiment of this invention. It will 
be apparent, however, that various variations and modifications may be made to the invention, with 
the attainment of some or all of the advantages of the invention. It is the object of the appended 
claims to cover these and such other variations and modifications as come within the true spirit and 
scope of the invention. 

What is claimed as new and desired to be secured by Letters Patent of the United States is: 
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